Y-Scaling Analysis of the Deuteron Within the Light-Front Dynamics Method 



O 

o 



oo 



O 

o 



X 



M.K. Gaidarov, M.V. Ivanov, and A.N. Antonov 

Institute of Nuclear Research and Nuclear Energy, 
Bulgarian Academy of Sciences, Sofia 1784, Bulgaria 

The concept of relativistic scaling is applied to describe the most recent data from inclusive 
electron-deuteron scattering at large momentum transfer. We calculate the asymptotic scaling 
function f(y) of the deuteron using its relationship with the nucleon momentum distribution. 
The latter is obtained in the framework of the relativistic light-front dynamics (LFD) method, 
in which the deuteron is described by six invariant functions fi (i=l,...,6) instead of two (S and 
D waves) in the nonrelativistic case. Comparison of the LFD asymptotic scaling function with 
other calculations using S and D waves corresponding to various nucleon-nucleon potentials, 
as well as with the Bethe-Salpeter result is made. It is shown that for \y\ > 400 MeV/c the 
differences between the LFD and the nonrelativistic scaling functions become larger. 



I. INTRODUCTION 



High-energy electron scattering from nuclei can provide important information on the wave function of nu- 
cleons in the nucleus. In particular, using simple assumptions about the reaction mechanism, scaling functions 
y—i ■ can be deduced that, if shown to scale (i.e., if they are independent of the momentum transfer), can provide 
' information about the momentum and energy distribution of the nucleons. Several theoretical studies have 
indicated that such measurements may provide direct access to studies of short-range nucleon-nucleon (NN) 
1 correlation effects. 

Since West's pioneer work Q , there has been a growth of interest in y-scaling analysis, both in its experimental 
and theoretical aspects. This is motivated by the importance of extracting nucleon momentum distributions 
from the experimental data. West showed that in the impulse approximation, if quasielastic scattering from 
■ a nucleon in the nucleus is a dominant reaction mechanism, a scaling function F(y) could be extracted from 
' the measured cross section which is related to the momentum distribution of the nucleons. In the simplest 
approximation the corresponding scaling variable y is the minimum momentum of the struck nucleon along the 
' direction of the virtual photon. In principle, the scaling function F(y, Q 2 ) depends on both y and momentum 
transfer Q 2 (-Q 2 is the square of the four- momentum transfer), but at sufficiently high Q 2 values the dependence 
on Q 2 should vanish yielding scaling. However, any attempt to extract the momentum distribution from the y- 
scaling in electron-nucleus scattering faces the problem of estimating both effects from the final-state interactions 
(FSI) of the struck nucleon with the rest of the nucleus and from the nucleon binding. Previous calculations 
[R-pf suggest that the contribution from the final-state interactions should vanish at sufficiently high Q 2 . The 
FSI lead to sizable scaling violation effects only at low values of the three- momentum transfer |q| = \fQ 2 + v 2 



(y is the energy loss) [9 10 . The most important dynamical effects, such as binding corrections, which represent 
the fact that for complex nuclei the final spectator A — 1 system can be left in all possible excited states 
including the continuum, have been treated in in terms of spectral functions. This problem has been solved 
in pd| , ^2[ by introducing a new scaling variable which gives direct, global and independent of A link between 
the experimental data and the longitudinal momentum components. 

Recently inclusive electron scattering has been studied at the Thomas Jefferson National Accelerator Facility 
(TJNAF) with 4.045 GeV incident beam energy from C, Fe and Au targets [||| to Q 2 w 7 (GeV/c) 2 . Data 
were also taken using liquid targets of hydrogen and deuterium The data presented in p3| , pd| represent 
a significant increase of the Q 2 range compared to previous SLAC measurement p5| , in which an approach to 
the scaling limit for heavy nuclei is suggested but for low values of \y\ < 0.3 GeV at momentum transfers up 
to 3 (GeV/c) 2 , and, moreover, a scaling behaviour is observed for the first time at very large negative y (y= 
-0.5 GeV/c). From theoretical point of view the extended region of y measured at TJNAF is of significant 
importance since this is a regime where the nucleon momentum distribution is expected to be dominated by 
short-range NN correlations. On the other hand, it is interesting to note that contributions from short-range 
FSI may also result in a scaling-like behavior due to the small Q 2 dependence of these effects, and that these 
contributions are also dominated by short-range correlations. Obviously, a complete understanding of this 
electron-nucleus scattering requires a relativistic approach to the quantities related to the y-scaling analysis for 
a detailed comparison with the experimental data. 

A relativistic y-scaling has been considered in |l6| ] by generalizing the nonrelativistic scaling function to the 
relativistic case. Realistic solutions of the spinor-spinor Bethe-Salpeter (BS) equation for the deuteron with 
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realistic interaction kernel were used for systematic investigation of the relativistic effects in inclusive quasielastic 
electron-deuteron scattering. The approach of y-scaling presented in |l6f| is fully covariant, with the deuteron 
being described by eight components, namely the 3 S^ + , 3 D+ + , 3 S 1 , 3 D 1 , 1 -P 1 + ~, 1 P± + , 3 Pi~, and 3 Pf + 
waves. The first two waves directly correspond to the S and D waves in the deuteron, with the waves with 
negative energy vanishing in the nonrelativistic limit. It has been demonstrated in ]l6| that, if the effects 
from the negative energy P-states are disregarded, the concept of covariant momentum distribution can be 
introduced. 

Recently a successful relativistic description of the nucleon momentum distribution in deuteron has been 
done |L7| within the light-front dynamics method |ll|[l{|. The most important result from the calculations in 
|jl7f is the possibility of the LFD method to describe simultaneously both deuteron charge form factors (that 
has been shown in J2(J) and the momentum distribution. It is shown in |2l],[l^] that after the projection of 
the Bcthe-Salpeter amplitude on the light front, the six components of the LFD deuteron wave function are 
expressed through integrals over the eight components of the deuteron Bethe-Salpeter amplitude. Provided the 
nucleon-nucleon interaction is the same, these approaches incorporate by different methods the same relativistic 
dynamics. The wave functions in LFD are the direct relativistic generalization of the nonrelativistic ones in 
the sense that they are still the probability amplitudes. Therefore they can be used in the relativistic nuclear 
physics (e.g. @). 

The aim of our work is using the nucleon momentum distribution n(k) obtained with the LFD method to 
calculate the deuteron scaling function. The result for the asymptotic function is compared with the recent 
TJNAF data measured at six values of Q 2 . In particular, the scaling behavior observed for very large negative 
y providing momenta higher than those corresponding to existing experimental data for n(k) may allow to 
distinguish the properties of the covariant LFD method from the potential approaches. The comparison with 
the BS result for the scaling function serves as a test for the consistency of both covariant approaches treating 
the deuteron relativistically in the case of y-scaling. 

The paper is organized as follows. In Section II the definition and the physical meaning of the scaling variable 
and the scaling function are briefly reviewed together with some basic relations between the nucleon spectral 
function, the scaling function and the momentum distribution. The results for the nucleon momentum distri- 
bution in deuteron obtained within the LFD method are given in Section III. The calculated LFD asymptotic 
scaling function of the deuteron is presented in Section IV, where a comparison with the Bethe-Salpeter result 
and with some nonrelativistic calculations is also done. The summary of the present work is given in Section V. 



II. BASIC RELATIONS IN THE ^-SCALING METHOD 



The scaling function is defined as the ratio of the measured cross section to the off-shell electron-nucleon cross 
section multiplied by a kinematic factor: 

where Z and N are the number of protons and neutrons in the target nucleus, respectively, a p and a n are the 
off-shell cross sections, and M is the mass of the proton. In analysing quasielastic scattering in terms of the 
y-scaling a new variable y = y(q, v) is introduced. Then the nuclear structure function which is determined 
using the spectral function P(k, E) as 

rE max {q,v) /•k ma j,(q,i',E) 

F(q,v) =2tt / dE kdkP(k,E), (2) 

can be expressed in terms of q and y rather than q and v (see Eq. (Q)). In Eq. (g) E = E rn i„ + E\_ x is the 
nucleon removal energy with E^_ 1 being the excitation energy of the final A — I nucleon system. 

The most commonly used scaling variable y is obtained p2[ ] starting from relativistic energy conservation, 
setting fc = y, = 1 and the excitation energy £^^=0, and is defined through the equation 

v + M A = (M 2 + q 2 + y 2 + lyq) 1 ' 2 + {M 2 A _ X + y 2 ) 1 / 2 , (3) 

where Ma is the mass of the target nucleus and Ma-x is the mass of the A — 1 nucleus. Therefore, y represents 
the longitudinal momentum of a nucleon having the minimum removal energy (E — E m i n , i.e. P41 =0). 
At high values of q a pure scaling regime is achieved, where k min « \y — [E — E min )\ and Eq. (0) becomes 

/>oo />oo 

F(q,y) -> f(y) =2tt / dE kdkP(k,E). (4) 

JEmin •> '\y-(E-E min )\ 
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In Eq. (Q) the particular behavior of P(k, E) at large k and E is used in order to extend the upper limits of 
integration to infinity Q . 

In the deuteron one always has E* A _ l =0 1 so that the spectral function is entirely determined by the nucleon 
momentum distribution n(k), i.e. P(k,E) — n(k)6(E — E m i n ), and, consequently, k m i n = \y\ for any value of 
q. The scaling function (^) reduces to the longitudinal momentum distribution 

f(y) = 2tt / kdkn(k). (5) 
J\y\ 



III. NUCLEON MOMENTUM DISTRIBUTION IN THE LIGHT-FRONT DYNAMICS METHOD 

The relativistic deuteron wave function (WF) on the light-front ^(fc, n) depends on two vector variables: i) the 
relative momentum k and ii) the unit vector n along lu which determines the position of the light-front surface. 
Due to this, the WF is determined by six invariant functions fi (i=l,...,6) instead of two (S- and Z?-waves) in 

the nonrelativistic case. Each one of these functions depends on two scalar variables k and z = cos{k 1 nV In 
LFD these six functions are calculated within the relativistic one-boson-exchange model. As shown in |18| , in 
the nonrelativistic limit the functions fz-G become negligible, /i^ do not depend on z and turn into S- and 
£>-waves (/i~its, f2~ — uo) and the wave function ^(fc, n) becomes the usual nonrelativistic wave function. One 
of the most important properties of the functions found in Jl8| is that for k > 2-^2.5 fm -1 the component 
/g (being related mainly to 7r-exchange) exceeds sufficiently the S- and D-w&ves. This fact is very important 
in the calculations of n(k) in deuteron as it will be shown below. 

The LFD calculations have shown (for more details, see Ref. |L7|) that, as expected, the most important 
contributions to the total n(fc) give terms related to the /i, and /g functions 

n{k)~n x {k) + n 2 (k) + n 5 (k). (6) 

The contributions of Hi, 112, nyi — ni + ^2 and rig are compared in Fig. [l| It can be seen that, while the 
functions fi and f% give a good description of the y-scaling data of n(k) for k < 2 fm _1 (like the S- and D-wave 
functions in the nonrelativistic case), it is impossible to explain the high- momentum components of n{k) at 
k > 2 fm _1 without the contribution of the function f$. We note that the deviation of the total n(k) from the 
sum ni2 = ni+n2 starts at k around 1.8 fm _1 . All this shows the important role of NN interactions which 
incorporate exchange of relativistic mesons in the case of the deuteron. 




1 2 3 4 5 



k [fm 1 ] 

FIG. 1. The nucleon momentum distribution in deuteron. The contributions of m, ri2, ni2 = m + 712 and 715 are 
presented. The y-scaling data are from Q|. The normalization is: J n(k)d 3 k = 1. 



The nucleon momentum distribution in the deuteron using S- and £>-wave functions ^s(k) and \E , £>(fc) 
corresponding to various NN potentials, such as the charge-dependent Bonn potential |23], the Argonne vig 
p4j| , the Nijmegen - I ,- II and - Reid 93 (|^] and Paris 1980 |Q i s a l so computed in (l7) by the expression: 

= ^s(k) + ^l(k)] = n s {k) + n D {k) (7) 
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with 



n(k)d 3 k = 1 . 



(8) 



In Fig. U the result for n(k) using the charge-dependent Bonn potential |^3) is given and compared with the 
y-scaling data. As can be seen, the D-component of n(k) is important but even its inclusion does not give a very 
good agreement with the data for k > 2 fm -1 . In the next Fig. ||we present the LFD result for n(k) compared 
with the calculations using the WF's corresponding to Nijmegen-I,-II, -Reid 93, Argonne vis and Paris 1980 
NN potentials and with the y-scaling data. We would like to note that: i) the results of the calculations using 
the NN potentials, such as Nijmegen-II, -Reid 93, Argonne Vis and Paris 1980 (shown in Fig. 0) are in better 
agreement with the y-scaling data than those using the charge-dependent Bonn potential (Fig. . This might 
be related to the fact that these potentials describe NN phase shifts up to larger energies (e.g. the Nijmegen-II 
potential gives reasonable pp phase shifts up to 1.2 GeV, while the charge-dependent one-boson exchange Bonn 
potential fits the phase-shift data below 350 MeV) ; ii) It can be seen from Fig. || that there are small differences 
between the curves corresponding to different NN potentials for k < 3 fm _1 (which give a good description of 
the y-scaling data and almost coincide with the LFD result) and larger ones for k > 3 fm _1 . Large differences 
take place, however, between all of them and the LFD result for k > 3 fin -1 . 




FIG. 2. The nucleon momentum distribution in deuteron (solid line) calculated using Eqs. (fj]) and (^) with S- and 
D-wave functions corresponding to the charge-dependent Bonn potential [^3|. S- and D-contributions are given by 
dashed and dotted line, respectively. The y-scaling data are from M. 




k [fm" 1 ] 

FIG. 3. The nucleon momentum distribution in deuteron calculated with the LFD (solid line) in comparison with the 
results presented by numerated arrows, as follows: 1- Argonne vig, 2-Nijmegen Reid 93, 3-Nijmegen I, 4-Nijmegen II, 
5-Paris 1980 NN potentials and with the y-scaling data M. 
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IV. RESULTS FOR THE ASYMPTOTIC SCALING FUNCTION OF DEUTERON 



In this Section we present the results for the asymptotic scaling function of deuteron which is calculated by 
using of the LFD nucleon momentum distribution given in Section III. 

The scaling function for deuteron calculated within the LFD method is shown in Fig. It is compared with 
the TJNAF experimental data |Q for all measured angles. The Q 2 values are given for Bjorken x = Q 2 /2Mv = 
1 and correspond to elastic scattering from a free nucleon. It is seen from Fig. |j that the relativistic LFD scaling 
function is in good agreement with the data in the whole region of negative y available. As known, the scaling 
breaks down for values of y > due to the dominance of other inelastic processes beyond the quasielastic 
scattering. Our LFD deuteron scaling function is also compared in Fig. || with the scaling function obtained 
within the BS formalism A small difference between the two results is observed for y < —400 MeV/c but, 
at the same time, the theoretical LFD scaling function is closer to the experimental data in the same region 
of y. The fact that both LFD and BS functions reveal similar behavior is a strong indication in favor of the 
consistency of the two relativistic covariant approaches in case of the y-scaling. 
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FIG. 4. The scaling function of deuteron. The experimental data for different values are from [Q. The solid and 
dashed curves represent the LFD calculations of this work and BS result of Ref. fl6| . 



In Fig. |5| the asymptotic relativistic LFD and BS scaling functions f(y) are compared with the nonrelativistic 
ones, calculated with some realistic interactions. We have already shown (see Figs. || and || of the present 
work) that the results using these potentials explain very well (with the exception of the charge-dependent 
Bonn potential) all the available data for n(k) up to k ~ 3 fm -1 exactly like the LFD method. However, it is 
concluded in that for k > 3 fm -1 the LFD results for n(k) deviate strongly from those of the calculations 
using NN potentials. Here we would like to emphasize the existence of the same discrepancies between the 
scaling functions observed from the comparison in Fig. 0. It is shown that for |y| > 400 MeV/c both LFD and 
BS curves start to deviate from the nonrelativistic scaling functions. The result for f(y) calculated using the 
Nijmegen-II NN potential is in better agreement with the experimental data than those using other potentials. 
It is in accordance with the result for n(k) shown in Fig. |3[ For instance, by a thorough comparison between 
the relativistic Bethe-Salpeter and the nonrelativistic scaling functions of deuteron it has been found in |l^] that 
the two functions start to sensibly differ also at \y\ > 400 MeV/c. Thus, the necessity to treat realistically the 
relativistic dynamics inside the deuteron in a way different from the potential approaches becomes apparent. 
In this sense, the results calculated for both momentum distribution and asymptotic scaling function confirm 
the abilities of the LFD method to describe with a good accuracy the experimental data measured at high 
momentum transfers. 
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FIG. 5. The asymptotic scaling function of deuteron. Line convention refering to calculations using LFD and potential 
approaches, as well as BS result ]lq| is given (see also the text). The experimental data [E7j are given by the full circles. 



V. CONCLUSIONS 

In the present paper inclusive electron-deuteron scattering data have been analyzed in terms of the y-scaling 
function within the light-front dynamics method. For this purpose, the nucleon momentum distribution in 
deuteron has been used in order to calculate the asymptotic scaling function. For the trivial case of deuteron, 
for which the structure function (Eq. ^) coincides with the longitudinal momentum distribution (Eq. ||) we have 
found a good agreement of the calculated scaling function with the experimental data. Thus, the concept of 
relativistic y-scaling can be introduced in the LFD relativistic description of inclusive quasielastic eD scattering, 
in the same way as it is done in the conventional nonrclativistic approach, i.e. by introducing a scaling function 
(which, in the scaling regime, is nothing but the nucleon longitudinal momentum distribution), and in terms of 
the same variable y. It has been pointed out that for \y\ > 400 MeV/c the differences between the LFD and 
the nonrelativistic scaling functions are very large. 

Exploring the light-front dynamics, we continue in this paper our analysis of important deuteron characteris- 
tics. The effective inclusion of the relativistic nucleon dynamics and of short-range NN correlations can be better 
seen when analyzing electron scattering at high momentum transfer from complex nuclei, for which a proper 
theoretical y-scaling analysis is still lacking. Although scaling violation effects due to final-state interactions 
sharply decrease with increasing momentum transfer, a consistent treatment of both FSI and nucleon binding 
must be made in order to perform a precise comparison with the new TJNAF data. Such an investigation is in 
progress. 
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